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ABSTRACT 

The presence of young stars in the immediate vicinity and strong tidal field of SgrA* remains unexplained. 
One currently popular idea for their origin posits that the stars were bused in by an Intermediate-Mass Black Hole 
(IMBH) which has inspiraled into the Galactic Center a few million years ago. Yu and Tremaine (2003) have 
argued that in this case some of the old stars in the SgrA* cusp would be ejected by hard gravitational collisions 
with the IMBH. Here we derive a general expression for the phase-space distribution of the ejected high-velocity 
stars, given the distribution function of the stars in the cusp. We compute it explicitly for the Peebles-Young 
distribution function of the cusp, and make a detailed model for the time-dependent ejection of stars during the 
IMBH inspiral. We find that (1) the stars are ejected in a burst lasting a few dynamical friction timescales; if the 
ejected stars are detected by Gaia they are likely to be produced by a single inspiral event, (2) if the inspiral is 
circular than in the beginning of the burst the velocity vectors of the ejected stars cluster around the inspiral plane, 
but rapidly isotropise as the burst proceeds, (3) if the inspiral is eccentric, then the stars are ejected in a broad jet 
roughly perpendicular to the Runge-Lenz vector of the IMBH orbit. In a typical cusp the orbit will precess with 
a period of ~ 10 5 years, and the rate of ejection into our part of the Galaxy (as defined by e.g. the Gaia visibility 
domain) will be modulated periodically. Gaia, together with the ground-based follow-up observations, will be 
able to clock many high-velocity stars back to their ejection from the Galactic Center, thus measuring some of the 
above phenomena. This would provide a clear signature of the IMBH inspiral in the past 10 — 20 Myr. 

Subject headings: black holes, massive stars, orbits 



1. INTRODUCTION 

The origin of young stars near SgrA* has been a major puzzle 
of the Galactic Center (GC) for over a decade (Phinney 1989, 
Sanders 1992, Morris 1993). One of the few promising ex- 
planations assumes that the stars were originally bound to an 
Intermediate-Mass Black Hole (IMBH) which has inspiraled 
from large distances (few parsecs) to the immediate vicinity 
(~ O.Olpc) of SgrA* (Hansen and Milosavljevic, 2003). In this 
scenario the IMBH itself would be a compact remnant of a su- 
perstar of several thousand solar masses. This superstar could 
be a product of stellar mergers at the center of dense young mas- 
sive cluster (Portegies Zwart et. al. 2004, Gurkan et. al. 2004) 
or could be produced in a fragmenting massive accretion disc 
(Levin and Beloborodov 2003, Levin 2003, Goodman and Tan 
2004). 

Regardless of its origin, the inpiraling IMBH will eject a 
number of stars from the GC and send them flying through (and 
out of) the Galaxy at high velocities (Yu and Tremaine 2003, 
see also Gualandris et. al. 2005). Brown et. al. 2005 have re- 
cently identified a high-velocity star in SDSS data and have ar- 
gued that it may have been ejected from the GC. The European 
satellite Gaia, scheduled to fly in 2012, will detect and measure 
proper motions of stars to magnitude 20; the Sun has magni- 
tude of 19.6 at lOkpc from Earth (Perryman et. al. 2001). The 
threshold for radial velocity measurement by Gaia is higher- 
magnitude 17; however the few high-proper-motion stars de- 
tected by Gaia could be observed by the ground-based tele- 
scopes and their radial velocities could be determined. Thus, 
many high proper-motion-velocity stars could potentially be de- 
tected, their spatial position in the Galaxy could be measured 
and their radial velocities could be determined from follow- 
up observations. With their 6-d phase-space information de- 



termined, their trajectories could be traced back to the GC and 
their times of ejection could be approximately determined; this 
point was recently and independently emphasized by Gnedin 
et. al. (2005). In this paper we derive the phase-space distri- 
bution of the ejected stars, and explain qualitatively how 6-d 
information could be used to constrain the IMBH orbital evolu- 
tion. In the next paper we will carry out the rigorous modeling 
of the signal expected by GAIA from an inspiraling IMBH. 

2. DISTRIBUTION OF EJECTED STARS 

The typical mass ratio of an IMBH and SgrA* ranges from 
10~ 2 to 10~ 4 . Therefore the gravitational interaction of a star 
with the IMBH which results in the star's ejection occurs very 
close to the IMBH and can be considered instantaneous, in a 
sense that the IMBH velocity is not affected by the central black 
hole (SgrA*) during the interaction. Moreover, the IMBH has 
an immediate effect on the distribution function of stars in the 
cusp only in it's immediate vicinity; therefore the stars interact- 
ing with the IMBH have initial velocities drawn from the distri- 
bution function of the cusp at the instantaneous location of the 
IMBH. The differential cross-section of instantaneous encoun- 
ters and the associated rates of stellar ejections were derived 
by Michel Henon in a series of seminal papers (1960a,b, 1969, 
see also Lin and Tremaine 1980). Our calculations are different 
from Henon's, since we care about the distribution of the veloc- 
ity vectors of ejected stars, and not just their energies. Also, we 
consider a single ejector (the IMBH) and do not average over 
the population of ejectors. The derivations below are of geo- 
metric nature; they are informed by but do not follow Henon's 
work. 

Consider an IMBH flying through the stellar cusp with ve- 
locity Vbh at radius r from SgrA*. It is reasonable to assume 
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that the cusp is composed of stars bound to the central black 
hole. Therefore, in velocity space, all the stars interacting with 
the IMBH are inside the sphere bounded by the escape velocity 
at r, 



2GM 



(1) 



see Fig. 1. In the equation above we have assumed that the 
potential is dominated by central black hole of mass M, which 
is an excellent approximation at the location of the puzzling 
young stars. 

First, let us calculate the rate of ejection of stars with velocity 
v > v e j, where 



v e j = q x v e 



(2) 



and q > 1. Here and further v, v esc , etc. stand for the magni- 
tudes of the corresponding velocity vectors. Figure 1 illustrates 
the velocity-space of the problem at hand. Let Vj n and v out be 
the velocity vector of a star before and after scattering by the 
IMBH, respectively. It is convenient to consider the problem in 
the frame of reference of the IMBH, and define 



Win = V fa -Vbh, 
WW = Vout-Vbh, 

W = I W in I = |v?out| • 



(3) 

(4) 



The dashed region of the velocity space marks the stars 
which can be scattered to velocities v out > v e j; analytically it 
is defined by two conditions: 

Vin ^ V esc j 

w > v e j-v bh . (5) 
The vector of the outgoing relative velocity w out must satisfy 

|vbh + W out | > V e j, (6) 

which can be translated into a constraint on the angle <fi out be- 
tween Vbh and w out : 



cos out > cos 00 = 



22 2 
2v b hW 



(7) 



The rate of ejection with velocities greater than v e j is given by 

#(V ut > V e j) = / d 3 V in f(v in )wA(w in ,(j)o), (8) 



where the integral is performed over the dashed volume in 
Fig. 1, f(v) is the velocity distribution function, and A(w m ,(f>o) 
is the cross-section for the star with initial relative velocity w;,, 
to be scattered into the velocity cone described by Eq. Q. This 
cross-section is evaluated in the Appendix: 



2x2 



A(w in , </>(>) = 7r(Gm/w ) 



sin 2 4>q 



(cOS0j n -COS0o) 2 ' 



(9) 



where m is the mass of the IMBH and </>;„ is the angle between 
Vbh and Win. 

Now, let us find the differential ejection rate dR(v oul )/dV , 
where dV is the incremental volume of the velocity space cen- 
tered on v out and dR is the rate of ejections of stars into this vol- 
ume. It is convenient to express dV as w 2 x dw x d£l oat where 
w = I Wout I = I Vout -v*bh I and c/£l ut is the increment of the solid 



angle around the vector w ou t- The differential ejection rate is 
then given by 

dR(v out )/dV = w / m b +Wn)dn J a( f« M . (10) 

J fl"out 

Here da(wi B , w ou t) / o^out is the differential crossection for grav- 
itational scattering from Wj n to w ou t, and c/f2i„ is the solid angle 
increment of the integrated vector variable Wj n . The differen- 
tial ejection rate is particularly easy to evaluate if the distribu- 
tion function /(vln) is axially symmetric about the vector Vbh- 
The ejection rate is then also axially symmetric. By averag- 
ing Eq. J over velocity vectors v oat which are connected by 
rotation about v b h one can show that 



dR(v oat )/dV = w / /(w, in ) sin 4> 



9A(W, <pm,(/) ut) 
1 S 1 ' 

dcos^out 



where 



COS 01 



v — vr. — w 

esc y bh ' 

2wv bh 



(ID 



(12) 



Here /(w, 4> in ) = /(v bh + w in ), and A(w, in , 0out) = A(w in , out ) 
as given by Eq. 0. 

For a general distribution function f(v), the integrals in 
Eqs. l|8) and Jl li must be integrated numerically. In the next 
section, we consider an important particular case where these 
integrals can be evaluated analytically. 

3. IMBH INSPIRALING THROUGH THE p <x R~ 15 CUSP 

A commonly considered cusp distribution function which is 
convenient for analytical work is the Peebles-Young distribu- 
tion function f(v) = /o; this corresponds to an isotropic cusp 
with the density profile p oc r" 1 5 (Young 1980). This profile is 
also very close to that of the observed cusp in the Galactic Cen- 
ter (Genzel et. al. 2003; see below), and hence in this section we 
will consider in some detail the distribution of the stars ejected 
by the IMBH inspiraling through the Peebles-Young cusp. 

The ejection rate is obtained by evaluating the integral in 
Eq. JSJ. We have found w and </>i n to be convenient integra- 
tion variables; it is easiest to integrate first over and then 
over w. One then finds the following for the ejection rate: 

3tt n(Gm) 2 . 

fl(Vout > V ej ) = — ^-3-^A<>,3), (13) 

^ v esc 

where p = v b h/v esc < 1, q = v e j/v eS c > 1, n is the local number 
density of stars and 



A(p.q) = 



a+2p-q) 



p 2 + q 2 1 
q 2 -l ~2 
2 2 r 

q -p 



(14) 



q-p 1+p 



2p 



2 21 

q -p 1 



q 2 -\ 



(l+pY-iq-pY fl+p 
-In 1 



6p(q 2 -l) 



q-p 



The differential ejection rate is then obtained by evaluating the 
integral in Eq. and is given by 

jrw -» s/, T , ?r/(Gm) 2 f 
dR(v oul )/dV = 5 {2C 



w 3 
1-C 2 



1 



C-C, C+l 



1 



(C-d) 2 (C+l) 2 



(15) 
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where 



W 
C 

Ci 



'bi, 



-2v out v bh cos9, 



(v out cos6»-v b h)/w, 



(16) 



2 2 2 

v - vr. — w 

esc ^bh ' 



2wv 



bh 



bh- 



and 9 is the angle between v out and v 

Let us put make a numerical estimate for Eq. Jl 3I > with the 
parameters relevant for the Galactic Center. We are interested 
in stars with luminosity no less than solar, since we want the 
star to be observable by Gaia at the distance of ~ 10 kpc. At 
the moment the observational constraint on the number of such 
stars in the central cusp is poor, as the faintest stars observ- 
able by Genzel's and Ghez' groups with the VLT and Keck 
have masses > 5M (R. Genzel and A. Ghez 2005, private 
communications). Among the stars which are observable, the 
mass function does not follow Salpeter's law, and is clearly top- 
heavy. The IMF favoring high masses is consistent both with 
stars born in a dense accretion disc (Artymowicz et. al. 1993, 
Levin 2003, Levin and Beloborodov 2003) or brought in as a 
part of a core of a dense young cluster (Gerhard 2001). To 
make further progress we need to parametrize the distribution 
of stars with masses ~ M© in the central cusp: 



n(r) = 8 x 10 6 pc~ 3 "o(0.1pc/r) 



1.5+5 



(17) 



The slope 6 was measured by Genzel et. al. (2003) and Shodel 
et. al. (2005) to be -0.2 and 0.2, respectively, with both of the 
analyses assuming constant mass-to-light ratio. Thus the cusp 
slope is close to the Peebles-Young profile, and in the fiducial 
example below we will assume 5 = 0. As explained above, the 
normalization factor no is not well known at present; no = 1 
corresponds to the case when the mass density at O.lpc quoted 
in Eq. (2) of Genzel et. al. (2003) is dominated by solar-mass 
stars. With this parametrization the ejection rate becomes 



tf(v ou t > v e j) ^ 60Myr l n Q x 



m/M 
0.001 



M 



1/2 



y 4 x 1O 6 M 

(18) 

Consider now an IMBH of mass m = 10 _3 M = 4000M Q which 
is inspiraling into SgrA* on a circular orbit. For a circular 
IMBH orbit, p = 1 /a/2 and the function A(l /y/2,q) is plotted 
as a function of q in Fig. 2. For the Peebles-Young profile the 
inspiral proceeds at constant rate, and the radius and the escape 
velocity of the orbit are given by 



A(p,q). 



r = r exp(-f /faction), 



= y/ 2GM/ r exp(f /2f frictioll ) 



where 



^friction — 



M 



1.5 



y/Gmpr l5 Q 



10 6 yr x (5/h () 



(19) 



(20) 



is the dynamical friction timescale, with Q being the Coulomb 
logarithm [cf. Eqs. (13) and (14) of Levin, Wu, and Thommes 
2005; our Q is their log A]. Let us now calculate, as a func- 
tion of time, the number of stars which are ejected during the 
inspiral and which will have velocities greater than v ga i axy = 



800km/sec once they climb through the bulge potential. We 
assume that the escape velocity from the bottom of the bulge is 
Vbuige = 300km/sec (see Yu and Tremaine 2003 for discussion). 
Therefore, after collision with the IMBH the star must have ve- 
locity greater than 



*esc * bulge galaxy 



(21) 



To evaluate the number of ejected solar-type stars as a function 
of time passed since the beginning of the inspiral, we integrate 
Eqs. Jl 81 . ( I19> , and \2\\ with respect to the time variable. The 
result is shown in Fig. 3, for a particular case of «o = 1 and Q = 5. 
The rate of ejections as a function of time is shown in Fig. 4. 
The solid part of the curve, i.e. the part with t < 3.25Myr, marks 
the first stage of the inspiral when the IMBH mass is smaller 
than the mass of the stellar cusp out to the radius of the IMBH 
orbit. Therefore, the Black-Hole inspiral does not significantly 
affect the stellar distribution function of the cusp, and our anal- 
ysis is reliable. For t > 3.25Myr, which corresponds to orbital 
radius < 0.015pc, the IMBH dominates the mass budget inside 
its orbit and depletes its neighborhood of stars on the timescale 
of ffriction- Therefore, within our analytical formalism we cannot 
reliably predict the rate of ejections since the original distribu- 
tion function of the cusp is strongly altered by the IMBH 1 . 
We model very roughly this IMBH-dominated regime by posit- 
ing that the overall density of stars and the rate of inspiral are 
declining exponentially on the timescale of ffriction, while as- 
suming that the overall form of the cusp distribution function 
remains the same. The dashed part of the curve shows the num- 
ber of ejections calculated within this rough model. To get a 
better intuition for the parameters of ejected stars, we have gen- 
erated mock data based on the ejection rates computed above. 
We present this data as a scatter plot in Fig. 5: on the x-axis we 
show the time of ejection and on the y-axis we show the veloc- 
ity of the ejected star after it escapes from the bulge. With Gaia 
one will be able to clock the high-velocity stars back to their 
ejections and thus measure both f e jection and 3-d velocity for 
many of the high- velocity stars within 10 kpc from the Earth 
whose proper-motion vector is pointing to the Galactic Cen- 
ter. Building a scatter plot like the one in Fig. 5 for the ejected 
high-velocity stars will provide a clear diagnostic for the IMBH 
inspiral as the origin of the ejected stars. 

On the theoretical side, we have a problem-more than half 
of the ejected stars are produced during the part of the inspi- 
ral which cannot be modeled analytically. Therefore, we have 
conceived, but not yet carried out, a series of numerical experi- 
ments aimed at addressing the part of the inspiral where IMBH 
mass dominates that of the stellar cusp. For now, we limit our- 
selves to two qualitative remarks: 

(a) The number of the stars ejected during a single inspiral is 
dependent on the poorly-known distribution of stars in the cusp 
but may be counted in hundreds. These stars are potential tar- 
gets for GAIA. 

(b) The stars are ejected in a burst the most intense part of which 
lasts 2 — 3 dynamical friction timescales. If no more than one 
IMBH per 10 Myr inspirals into SgrA*, than the high-velocity 
stars seen by GAIA should be all generated by a single inspiral 
event. 

How anisotropic are the ejected stars? First we must note that 
the velocity vector v star of a star after it escapes the potential of 



'Yu and Tremaine (2003) have estimated the rate of ejections in this regime by appealing to the loss-cone formalism. They have assumed that the BH binary 
interacts with the stars on plunging orbits and that the star with pericenter radius similar to that of the binary gets ejected almost immediately. The validity of both 
assumptions is unclear for the small (10~ 3 ) mass ratios considered here. 
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the central black hole is different from v ou t> which is the stel- 
lar velocity right after gravitational collision with the IMBH. 
Generally the velocity of the escaped star 



Vstar — V s tar('" ; v oul) 



(22) 



can be found by solving for the Keplerian hyperbolic trajectory 
of the star; here r is the position of the IMBH during the grav- 
itational collision. The mathematical expressions for Keplerian 
hyperbolae can be found in a standard textbook on Celestial 
Mechanics, e.g. Danby (1988). 

A useful way to quantify the anisotropy of ejections is to 
compute the dispersion of angles x(^star) between the velocity 
vectors of the ejected stars and the inspiral plane. This disper- 
sion cr can be calculated by evaluating the following integral 
over a single period of the IMBH orbit: 



a 2 = - 



1 



N ei (T) 



dt I d 3 v. 



dR{v bh (t),v 0Ut ) 
dV 



XX [v star (fU),Vout)], 

(23) 



where t is the time, T is the period of the orbit, dR / dV is given 
by Eq. (llOi [and in appropriate cases by Eq. il It or dl5H . and 
N e j(T) is the number of stars ejected during an orbital period: 



dt 



d 3 v c 



<//?(Vbh(0,Vout) 

dV 



(24) 



In Fig. 6 we plot a vs time for the circular IMBH inspiral, with 
the thick horizontal line corresponding to isotropic distribution 
of the velocity vectors. At early times the stars ejected with 
velocities > 800km/sec will be clustering towards the inspiral 
plane, while at later times the velocity vectors of the ejected 
stars have no strong anisotropy. Unfortunately, at early times 
the ejection rate is low, and thus in the example we are consider- 
ing only the first 10 or so ejections will show strong anisotropy. 
GAIA will see at most about 1/4 of the ejected stars; we think 
therefore that detection of spatial anisotropy for circular inspi- 
ral is unlikely. 

Eccentric inspiral. There is a stronger chance to detect spa- 
tial anisotropy of the ejected stars if the IMBH inspirals on 
a strongly eccentric orbit. Most high-velocity stars will be 
ejected when IMBH passes the pericenter. These stars will then 
form a broad jet which will have the general direction of the 
IMBH velocity vector at the pericenter, i.e. perpendicular to the 
Runge-Lenz vector of the IMBH orbit. The elliptical orbit of 
the IMBH will precess in the potential of the cusp, and the jet 
will rotate with the orbital ellipse. Therefore, the rate of stars 
ejected into the region of the Galaxy observable by GAIA (the 
GAIA sphere) will be modulated with the period of the IMBH 
precession. It is this periodic modulation that is potentially ob- 
servable, and that would be a tell-tale of the eccentric IMBH 
inspiral. 

For the Peebles- Young cusp the period of precession depends 
only on the IMBH eccentricity and does not depend on its semi- 
major axis. Thus during the first stage of the inspiral, when the 
cusp is not affected by the IMBH, the period of modulation will 
be constant and equal to 



'precession ' 



1.3 



M 



4 x 1O 6 M 



n x 10 5 yr, 



(25) 



see Munyaneza et. al. (1999). Here we have taken the IMBH 
eccentricity 2 e = 0.7. The Galactic Center cusp has spatial dis- 
tribution similar to the Peebles-Young profile, and therefore the 
period of precession will be nearly constant during the first 
stage of the inspiral. During the second stage, however, when 
the IMBH depletes the stellar density of the cusp, the preces- 
sion period becomes longer. Detection of this period change 
would clearly provide a strong probe of the cusp structure. In 
Fig. 7 we show the rate of stars injected into the GAIA visibility 
sphere as a function of the injection time (all computed for the 
fiducial cusp parameters described above). The star is defined 
to be in the GAIA visibility sphere if at the current time, taken 
to be lOMyr, the star is within lOkpc from the Earth. The orbit 
of the IMBH is taken to lie in the Galactic plane. We clearly 
see in Fig. 7 both the periodicity in the injection rate and the 
change in this periodicity when the IMBH begins to empty its 
local region of the cusp. 

4. CONCLUSIONS. 

In this paper we have derived the analytic expression for the 
time-dependent phase-space distribution of stars ejected from 
the Galactic Center as a result of the IMBH inspiral. We have 
discussed how Gaia, together with the follow-up ground-based 
radial velocity measurements, may be used to track the history 
of ejections from the Galactic Center; this point was recently 
and independently made by Gnedin et. al. 2005. We find that 
a single IMBH inspiral produces a burst in the rate of ejection 
which has a duration a few dynamical friction timescales. In the 
very beginning of the burst, the velocities of the ejected stars 
cluster around the inspiral plane, however this anisotropy will 
be hard to detect and most of the stars will be ejected isotropi- 
cally. 

If the IMBH orbit is eccentric, then the ejected stars form a 
broad "jet" roughly aligned with IMBH velocity in the pericen- 
ter. Because of the cusp potential, the IMBH orbit will precess 
with the period of ~ 10 5 years, and the rate of stellar ejections 
into the GAIA visibility zone will be strongly periodically mod- 
ulated. When the IMBH begins to modify stellar distribution in 
its neighborhood (i.e., when its mass becomes comparable to 
the mass of stars inside its orbit), the period of precession will 
increase, as shown in Fig. 7. 

Therefore, the burst in the ejection rate and potential modula- 
tion of this rate during the burst, if detected by GAIA, would be 
a clear signature that an IMBH has inspiraled into the Galactic 
Center in the past 10-20 Myr. 
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APPENDIX: EVALUATION OF THE CROSSECTION A 

In this Appendix we evaluate the crossection for a star with 
the initial velocity w m relative to the IMBH to scatter into the 
cone given by Eq. Q 



COS( 



> COS 00 : 



v 2 - 



"bh 



2v bh w 



(26) 



2 Previous calculations (e.g., in the Appendices of Gould and Quillen 2003 and Levin et. al. 2005) argue that the eccentricity of the IMBH inspiraling through the 
Peebles-Young cusp remains constant with time. However, these calculations neglect secular torques exerted on a precessing eccentric orbit; therefore the question 
of eccentricity evolution of an inspiraling IMBH remains open. 
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Here cf> oul is the angle between vv out and the velocity vector of 
the IMBH. 

Let P be the plane which is perpendicular to w m and which 
contains the IMBH. Let b be the impact parameter vector, i.e. a 
vector directed from the IMBH to the projection of the star onto 
P before the star is scattered by the IMBH. After the scattering 
is complete, the stellar velocity relative to the IMBH is given 

by 



Wont = Wi n cos 9 - b sin 9(w/b) , 
where 9 is the angle of deflection given by 

Gm 



tan(0/2) = 



w 2 b 



(27) 



(28) 



From Eq. (1271 we see that 



cos <f> out = cos 0i„ cos 9 + sin (f>- m sin 9 cos <f>b , 



(29) 



where </>/, is the polar angle in the fo-plane. 

Equations d29t and (1281 together with Eq. i26\ define a fi- 
nite domain in the b-plane. It is straightforward but somewhat 
tedious to find the boundary of this domain; this calculation is 
not shown here. The area of the domain is the crossection we 
seek and is given by 



Mwin, 4>o) = n(Gm/w 2 ) 2 



(COS <pi n - COS <po) 2 ' 

which is the Equation l|9} in the text. 



(30) 



REFERENCES 



Artymowicz, P., Lin, D. N. C, & Wampler, E. J. 1993, ApJ, 409, 592 
Brown, W. B., et. al. 2005, ApJ, 622, L33 

Danby, J. M. A. 1988, Fundamentals of Celestial Mechanics, Richmond, Va., 
USA, Willmann-Bell, 2nd ed. 

Gerhard, O. 2001, ApJ, 54 6, L39 

Gnedin, O. Y., et. al. 2005, astro-ph/0506739 
Goodman, J., & Tan, J. C. 2004, ApJ, 608, 108 

Gualandris, A., Portegies Zwart, S., & Sipior, M. 2005, accepted to MNRAS, 

astro-ph/0507365 
Gurkan, M. A., hreitag, M., & Rasio, F. A. 2004, ApJ, 604, 632 
Hansen, B. M. S., & Milosavljevic, M. 2003, ApJ, 593, L77 
Henon, M. 1960a, AnAp, 23, 467 
Henon, M. 1960b, AnAp, 23, 668 
Henon, M. 1969, A&A, 2, 151 
Levin, Y. 2003, astro-ph/0307084 



Levin, Y, & Beloborodov, A. M. 2003, ApJ, 590, L33 

Levin. Y. Wu. A S. P., Thommes, E. W. 2005, submitted to ApJ, 

astro-ph/0502143 
Lin, D. N. L'., & Tremaine, S. 1980, ApJ, 264, 364 
Morris, M. 1993, ApJ, 408, 496 

Munyaneza, F, Tsiklauri, D., Viollier, R. D. 1999, ApJ, 526, 744 

Phinney, E. S. 1989, in IAU Symp. 136, The Center of the Galaxy, ed. M. Morris 

(Dordrecht: Kluwer), 543 
Perryman, M. A. C, et. al. 2001, A&A, 369, 339 
Portegies Zwart, S. F, et. al. 2004, Nature, 428, 724 
Sanders, R. H. 1992, Nature, 359, 131 
Shodel, R., et. al. 2005, in preparation 
Young, P. 1980, ApJ, 242, 1232 
Yu, Q., & Tremaine, S. 2003, ApJ, 599, 1 129 



6 




Fig. l.— Velocity-space diagram of gravitational scattering of stars by the IMBH. The initial stellar velocities vl n lie inside the sphere 
of radius v esc . The stars with velocities in the shadowed region of the figure can be scattered by the IMBH to velocities greater than 
v e j. Other symbols are explained in the text. 
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Fig. 2.— The function A(p,q) computed for a circular IMBH orbit, p = 1 j\[2 
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Fig. 3.— The number of ejected stars as a function of the ejection time. The solid curve represents the region where analytical 
calculations are appropriate, the dashed curve represents the region where our model is a rough approximation. 
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Fig. 4.— The rate of ejections as a function of the ejection time. 
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Fig. 5.— Scatterplot of velocities and ejection times for the ejected stars 
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Fig. 6.— Dispersion of angles between the ejection velocity vectors and the inspiral plane plotted as a function of ejection time 




Fig. 7. — Eccentric inspiral: the rate of ejections of stars into GAIA visibility domain as a function of the ejection time. Here the 
IMBH orbital eccentricity is e = 0.7. 



